Comparative analysis of genes expression involved in type II toxin-antitoxin system in Staphylococcus aureus following persister cell formation.

BACKGROUND: The formation of persister cells is the main reason for persistent infections. They are associated with antibiotic treatment failure and subsequently chronic infection. The study aimed to assess the expression of type II toxin/antitoxin (TA) system genes in persister cells of Staphylococcus aureus in the presence of the following antibiotics vancomycin, ciprofloxacin, and gentamicin in exponential and stationary phases. METHODS AND RESULTS: The colony count was used to evaluate the effect of different types of antibiotics on S. aureus persister cell formation during exponential and stationary phases. Moreover, the expression level of TA systems and clpP genes in the persister population in exponential and stationary phases were measured by quantitative reverse transcriptase real-time PCR (qRT-PCR). The results of the study showed the presence of persister phenotype of S. aureus strains in the attendance of bactericidal antibiotics in comparison to the control group during the exponential and stationary phases. Moreover, qRT-PCR resulted in the fact that the role of TA systems involved in the persister cell formation depends on the bacterial growth phase and the type of strain and antibiotic. CONCLUSIONS: In total, the present study provides some data on the persister cell formation and the possible role of TA system genes in this process.

PCR Development for Analysis of Some Type II Toxin-Antitoxin Systems, relJK, mazEF3, and vapBC3 Genes, in Mycobacterium tuberculosis and Mycobacterium bovis.

Toxin-Antitoxin (TA) systems are some small genetic modules in bacteria that play significant roles in resistance and tolerance development to antibiotics. Whole genome sequencing (WGS) is an effective method to analyze TA systems in pathogenic Mycobacteria. However, this study aimed to use a simple and inexpensive PCR-Sequencing approach to investigate the type II TA system. Using data from the WGS of Mycobacterium tuberculosis (M. tuberculosis) strain H37Rv and Mycobacterium bovis (M. bovis) strain BCG, primers specific to the relJK, mazEF3, and vapBC3 gene families were designed by Primer3 software. Following that, a total of 90 isolates were examined using the newly developed PCR assay, consisting of 64 M. tuberculosis and 26 M. bovis isolates, encompassing both 45 rifampin-sensitive and 45 rifampin-resistant strains. Finally, 28 isolates (including 14 rifampin-resistant isolates) were sent for sequencing, and their sequences were aligned and compared to the mentioned reference sequences. The amplicons size of mazEF3, relJK, and vapBC3 genes were 825, 875, and 934 bp, respectively. Furthermore, all tested isolates showed the specific amplicons for these TA families. To evaluate the specificity of the primers, PCR was performed on S. aureus and E.coli isolates. None of the examined samples had the desired amplicons. Therefore, the primers had acceptable specificity. The results indicated that the developed PCR-Sequencing approach can be used to effectively investigate certain types of TA systems. Considering high costs of WGS and difficulty in interpreting its results, such a simple and inexpensive method is beneficial in the evaluation of TA systems in Mycobacteria.

Arbitrium communication controls phage lysogeny through non-lethal modulation of a host toxin-antitoxin defence system.

Temperate Bacillus phages often utilize arbitrium communication to control lysis/lysogeny decisions, but the mechanisms by which this control is exerted remains largely unknown. Here we find that the arbitrium system of Bacillus subtilis phage ϕ3T modulates the host-encoded MazEF toxin-antitoxin system to this aim. Upon infection, the MazF ribonuclease is activated by three phage genes. At low arbitrium signal concentrations, MazF is inactivated by two phage-encoded MazE homologues: the arbitrium-controlled AimX and the later-expressed YosL proteins. At high signal, MazF remains active, promoting lysogeny without harming the bacterial host. MazF cleavage sites are enriched on transcripts of phage lytic genes but absent from the phage repressor in ϕ3T and other Spß-like phages. Combined with low activation levels of MazF during infections, this pattern explains the phage-specific effect. Our results show how a bacterial toxin-antitoxin system has been co-opted by a phage to control lysis/lysogeny decisions without compromising host viability.

Cassette recombination dynamics within chromosomal integrons are regulated by toxin-antitoxin systems.

Integrons are adaptive bacterial devices that rearrange promoter-less gene cassettes into variable ordered arrays under stress conditions, thereby sampling combinatorial phenotypic diversity. Chromosomal integrons often carry hundreds of silent gene cassettes, with integrase-mediated recombination leading to rampant DNA excision and integration, posing a potential threat to genome integrity. How this activity is regulated and controlled, particularly through selective pressures, to maintain such large cassette arrays is unknown. Here, we show a key role of promoter-containing toxin-antitoxin (TA) cassettes as systems that kill the cell when the overall cassette excision rate is too high. These results highlight the importance of TA cassettes regulating the cassette recombination dynamics and provide insight into the evolution and success of integrons in bacterial genomes.

The Chlamydia-related Waddlia chondrophila encodes functional type II toxin-antitoxin systems.

Bacterial toxin-antitoxin (TA) systems are widespread in chromosomes and plasmids of free-living microorganisms, but only a few have been identified in obligate intracellular species. We found seven putative type II TA modules in Waddlia chondrophila, a Chlamydia-related species that is able to infect a very broad series of eukaryotic hosts, ranging from protists to mammalian cells. The RNA levels of Waddlia TA systems are significantly upregulated by iron starvation and novobiocin, but they are not affected by antibiotics such as ß-lactams and glycopeptides, which suggests different mechanisms underlying stress responses. Five of the identified TA modules, including HigBA1 and MazEF1, encoded on the Waddlia cryptic plasmid, proved to be functional when expressed in a heterologous host. TA systems have been associated with the maintenance of mobile genetic elements, bacterial defense against bacteriophages, and persistence upon exposure to adverse conditions. As their RNA levels are upregulated upon exposure to adverse conditions, Waddlia TA modules may be involved in survival to stress. Moreover, as Waddlia can infect a wide range of hosts including free-living amoebae, TA modules could also represent an innate immunity system to fight against bacteriophages and other microorganisms with which Waddlia has to share its replicative niche.IMPORTANCEThe response to adverse conditions, such as exposure to antibiotics, nutrient starvation and competition with other microorganisms, is essential for the survival of a bacterial population. TA systems are modules composed of two elements, a toxic protein and an antitoxin (protein or RNA) that counteracts the toxin. Although many aspects of TA biological functions still await to be elucidated, TAs have often been implicated in bacterial response to stress, including the response to nutrient starvation, antibiotic treatment and bacteriophage infection. TAs are ubiquitous in free-living bacteria but rare in obligate intracellular species such as chlamydiae. We identified functional TA systems in Waddlia chondrophila, a chlamydial species with a strikingly broad host range compared to other chlamydiae. Our work contributes to understand how obligate intracellular bacteria react to adverse conditions that might arise from competition with other viruses/bacteria for the same replicative niche and would threaten their ability to replicate.

dbAPIS: a database of anti-prokaryotic immune system genes.

Anti-prokaryotic immune system (APIS) proteins, typically encoded by phages, prophages, and plasmids, inhibit prokaryotic immune systems (e.g. restriction modification, toxin-antitoxin, CRISPR-Cas). A growing number of APIS genes have been characterized and dispersed in the literature. Here we developed dbAPIS (https://bcb.unl.edu/dbAPIS), as the first literature curated data repository for experimentally verified APIS genes and their associated protein families. The key features of dbAPIS include: (i) experimentally verified APIS genes with their protein sequences, functional annotation, PDB or AlphaFold predicted structures, genomic context, sequence and structural homologs from different microbiome/virome databases; (ii) classification of APIS proteins into sequence-based families and construction of hidden Markov models (HMMs); (iii) user-friendly web interface for data browsing by the inhibited immune system types or by the hosts, and functions for searching and batch downloading of pre-computed data; (iv) Inclusion of all types of APIS proteins (except for anti-CRISPRs) that inhibit a variety of prokaryotic defense systems (e.g. RM, TA, CBASS, Thoeris, Gabija). The current release of dbAPIS contains 41 verified APIS proteins and ∼4400 sequence homologs of 92 families and 38 clans. dbAPIS will facilitate the discovery of novel anti-defense genes and genomic islands in phages, by providing a user-friendly data repository and a web resource for an easy homology search against known APIS proteins.

TADB 3.0: an updated database of bacterial toxin-antitoxin loci and associated mobile genetic elements.

TADB 3.0 (https://bioinfo-mml.sjtu.edu.cn/TADB3/) is an updated database that provides comprehensive information on bacterial types I to VIII toxin-antitoxin (TA) loci. Compared with the previous version, three major improvements are introduced: First, with the aid of text mining and manual curation, it records the details of 536 TA loci with experimental support, including 102, 403, 8, 14, 1, 1, 3 and 4 TA loci of types I to VIII, respectively; Second, by leveraging the upgraded TA prediction tool TAfinder 2.0 with a stringent strategy, TADB 3.0 collects 211 697 putative types I to VIII TA loci predicted in 34 789 completely sequenced prokaryotic genomes, providing researchers with a large-scale dataset for further follow-up analysis and characterization; Third, based on their genomic locations, relationships of 69 019 TA loci and 60 898 mobile genetic elements (MGEs) are visualized by interactive networks accessible through the user-friendly web page. With the recent updates, TADB 3.0 may provide improved in silico support for comprehending the biological roles of TA pairs in prokaryotes and their functional associations with MGEs.

Preventing toxicity in toxin-antitoxin systems: An overview of regulatory mechanisms.

Toxin-antitoxin systems (TAs) are generally two-component genetic modules present in almost every prokaryotic genome. The production of the free and active toxin is able to disrupt key cellular processes leading to the growth inhibition or death of its host organism in absence of its cognate antitoxin. The functions attributed to TAs rely on this lethal phenotype ranging from mobile genetic elements stabilization to phage defense. Their abundance in prokaryotic genomes as well as their lethal potential make them attractive targets for new antibacterial strategies. The hijacking of TAs requires a deep understanding of their regulation to be able to design such approach. In this review, we summarize the accumulated knowledge on how bacteria cope with these toxic genes in their genome. The characterized TAs can be grouped based on the way they prevent toxicity. Some systems rely on a tight control of the expression to prevent the production of the toxin while others control the activity of the toxin at the post-translational level.

Unraveling the evolutionary dynamics of toxin-antitoxin systems in diverse genetic lineages of Escherichia coli including the high-risk clonal complexes.

IMPORTANCE: Large-scale genomic studies of E. coli provide an invaluable opportunity to understand how genomic fine-tuning contributes to the transition of bacterial lifestyle from being commensals to mutualists or pathogens. Within this context, through machine learning-based studies, it appears that TA systems play an important role in the classification of high-risk clonal lineages and could be attributed to their epidemiological success. Due to these profound indications and assumptions, we attempted to provide unique insights into the ordered world of TA systems at the population level by investigating the diversity and evolutionary patterns of TA genes across 19 different STs of E. coli. Further in-depth analysis of ST-specific TA structures and associated genetic coordinates holds the potential to elucidate the functional implications of TA systems in bacterial cell survival and persistence, by and large.

Structural basis for kinase inhibition in the tripartite E. coli HipBST toxin-antitoxin system.

Many bacteria encode multiple toxin-antitoxin (TA) systems targeting separate, but closely related, cellular functions. The toxin of the Escherichia coli hipBA system, HipA, is a kinase that inhibits translation via phosphorylation of glutamyl-tRNA synthetase. Enteropathogenic E. coli O127:H6 encodes the hipBA-like, tripartite TA system; hipBST, in which the HipT toxin specifically targets the tryptophanyl-tRNA synthetase, TrpS. Notably, in the tripartite system, the function as antitoxin has been taken over by the third protein, HipS, but the molecular details of how activity of HipT is inhibited remain poorly understood. Here, we show that HipBST is structurally different from E. coli HipBA and that the unique HipS protein, which is homologous to the N-terminal subdomain of HipA, inhibits the kinase through insertion of a conserved Trp residue into the active site. We also show how auto-phosphorylation at two conserved sites in the kinase toxin serve different roles and affect the ability of HipS to neutralize HipT. Finally, solution structural studies show how phosphorylation affects overall TA complex flexibility.

Effect of ZnO nanoparticles on biofilm formation and gene expression of the toxin-antitoxin system in clinical isolates of Pseudomonas aeruginosa.

BACKGROUND: Biofilm formation by Pseudomonas aeruginosa (P. aeruginosa) is known to be characteristic of this organism. This bacterium is considered one of the most life-threatening bacteria and has been identified as a priority pathogen for research by WHO. Biofilm-producing P. aeruginosa is a concern in many parts of the world due to antibiotic resistance. Alginate also plays an important role in the biofilm formation of P. aeruginosa as well as the emergence of antibiotic resistance in biofilms. In addition, the systems of toxin-antitoxin( TA) play an important role in biofilm formation. Metal nanoparticle(NP) such as zinc oxide (ZnO) also have extensive biological properties, especially anti-biofilm properties. Therefore, this study was conducted in relation to the importance of zinc oxide nanoparticles (ZnO NPs) in biofilm formation and also the correlation of gene expression of TA systems in clinical isolates of P. aeruginosa. METHODS: A total of 52 P. aeruginosa isolates were collected from burns (n = 15), UTI (n = 31), and trachea (n = 6) in hospitals in Ilam between May 2020 and October 2020. Biofilm formation was assessed using a microtiter plate assay. MIC and sub-MIC concentrations of ZnO NPs (10-30 nm with purity greater than 99.8%) in P. aeruginosa were determined. Subsequently, biofilm formation was investigated using sub-MIC concentrations of ZnO NPs. Finally, total RNA was extracted and RT- qPCR was used to determine the expression levels of genes of mazEF, mqsRA, and higBA of TA systems. RESULTS: Six isolates of P. aeruginosa were found to form strong biofilms. The results showed that ZnO NPs were able to inhibit biofilm formation. In our experiments, we found that the sub-MIC concentration of ZnO NPs increased the gene expression of antitoxins mazE and mqsA and toxin higB of TA systems treated with ZnO NPs. CONCLUSIONS: In the present study, ZnO NPs were shown to effectively inhibit biofilm formation in P. aeruginosa. Our results support the relationship between TA systems and ZnO NPs in biofilm formation in P. aeruginosa. Importantly, the expression of antitoxins mazE and mqsA was high after treatment with ZnO NPs, but not that of antitoxin higA.

Functional diversification despite structural congruence in the HipBST toxin-antitoxin system of Legionella pneumophila.

IMPORTANCE: Toxin-antitoxin (TA) systems are parasitic genetic elements found in almost all bacterial genomes. They are exchanged horizontally between cells and are typically poorly conserved across closely related strains and species. Here, we report the characterization of a tripartite TA system in the bacterial pathogen Legionella pneumophila that is highly conserved across Legionella species genomes. This system (denoted HipBSTLp) is a distant homolog of the recently discovered split-HipA system in Escherichia coli (HipBSTEc). We present bioinformatic, molecular, and structural analyses of the divergence between these two systems and the functionality of this newly described TA system family. Furthermore, we provide evidence to refute previous claims that the toxin in this system (HipTLp) possesses bifunctionality as an L. pneumophila virulence protein. Overall, this work expands our understanding of the split-HipA system architecture and illustrates the potential for undiscovered biology in these abundant genetic elements.

A novel double-ribonuclease toxin-antitoxin system linked to the stress response and survival of Acidovorax citrulli.

IMPORTANCE: Bacterial fruit blotch (BFB), which is caused by the seed-borne bacterium Acidovorax citrulli, is a devastating disease affecting cucurbit crops throughout the world. Although seed fermentation and treatment with disinfectants can provide effective management of BFB, they cannot completely guarantee pathogen-free seedstock, which suggests that A. citrulli is a highly stress-resistant pathogen. Toxin-antitoxin (TA) systems are common among a diverse range of bacteria and have been reported to play a role in bacterial stress response. However, there is currently much debate about the relationship between TA systems and stress response in bacteria. The current study characterized a novel TA system (Aave_1720-Aave_1719) from A. citrulli that affects both biofilm formation and survival in response to sodium hypochlorite stress. The mechanism of neutralization differed from typical TA systems as two separate mechanisms were associated with the antitoxin, which exhibited characteristics of both type II and type V TA systems. The Aave_1720-Aave_1719 system described here also constitutes the first known report of a double-ribonuclease TA system in bacteria, which expands our understanding of the range of regulatory mechanisms utilized by bacterial TA systems, providing new insight into the survival of A. citrulli in response to stress.

Toxin-Antitoxin system of Mycobacterium tuberculosis: Roles beyond stress sensor and growth regulator.

The advent of effective drug regimen and BCG vaccine has significantly decreased the rate of morbidity and mortality of TB. However, lengthy treatment and slower recovery rate, as well as reactivation of the disease with the emergence of multi-drug, extensively-drug, and totally-drug resistance strains, pose a serious concern. The complexities associated are due to the highly evolved and complex nature of the bacterium itself. One of the unique features of Mycobacterium tuberculosis [M.tb] is that it has undergone reductive evolution while maintaining and amplified a few gene families. One of the critical gene family involved in the virulence and pathogenesis is the Toxin-Antitoxin system. These families are believed to harbor virulence signature and are strongly associated with various stress adaptations and pathogenesis. The M.tb TA systems are linked with growth regulation machinery during various environmental stresses. The genes of TA systems are differentially expressed in the host during an active infection, oxidative stress, low pH stress, and starvation, which essentially indicate their role beyond growth regulators. Here in this review, we have discussed different roles of TA gene families in various stresses and their prospective role at the host-pathogen interface, which could be exploited to understand the M.tb associated pathomechanisms better and further designing the new strategies against the pathogen.

abkBA toxin-antitoxin system may act as antipersister modules in Acinetobacter baumannii clinical isolates.

Aim: Persistence cells comprise a subpopulation of bacteria that is resistant to treatment. In this study, the role of the toxin-antitoxin (TA) system in the formation of persistence cells of Acinetobacter baumannii isolates was investigated. Methods: After confirming all isolates, TA systems abkBA, mqsRA and higBA were identified. Persister cells were confirmed using the standard method. Real-time PCR was used to compare the expression of TA systems in isolates in persistence and normal states. Results: The abkAB system was present in all isolates; 4% of isolates formed persister cells. The expression level of the abkB gene in persistent isolates showed a sevenfold increase compared with nonpersistent isolates. Conclusion: The abkBA system is proposed as an antipersistence target in A. baumannii isolates.

Minimized antibiotic-free plasmid vector for gene therapy utilizing a new toxin-antitoxin system.

Approaches to improve plasmid-mediated transgene expression are needed for gene therapy and genetic immunization applications. The backbone sequences needed for the production of plasmids in bacterial hosts and the use of antibiotic resistance genes as selection markers represent biological safety risks. Here, we report the development of an antibiotic-free expression plasmid vector with a minimized backbone utilizing a new toxin-antitoxin (TA) system. The Rs_0636/Rs_0637 TA pair was derived from the coral-associated bacterium Roseivirga sp. The toxin gene is integrated into the chromosome of Escherichia coli host cells, and a recombinant mammalian expression plasmid is constructed by replacing the antibiotic resistance gene with the antitoxin gene Rs_0637 (here named Tiniplasmid). The Tiniplasmid system affords high selection efficiency (∼80%) for target gene insertion into the plasmid and has high plasmid stability in E. coli (at least 9 days) in antibiotic-free conditions. Furthermore, with the aim of reducing the size of the backbone sequence, we found that the antitoxin gene can be reduced to 153 bp without a significant reduction in selection efficiency. To develop its applications in gene therapy and DNA vaccines, the biosafety and efficiency of the Tiniplasmid-based eukaryotic gene delivery and expression were further evaluated in CHO-K1 cells. The results showed that Rs_0636/Rs_0637 has no cell toxicity and that the Tiniplasmid vector has a higher gene expression efficiency than the commercial vectors pCpGfree and pSTD in the eukaryotic cells. Altogether, the results demonstrate the potential of the Rs_0636/Rs_0637-based antibiotic-free plasmid vector for the development and production of safe and efficacious DNA vaccines.

Bacterial Toxin-Antitoxin Systems' Cross-Interactions-Implications for Practical Use in Medicine and Biotechnology.

Toxin-antitoxin (TA) systems are widely present in bacterial genomes. They consist of stable toxins and unstable antitoxins that are classified into distinct groups based on their structure and biological activity. TA systems are mostly related to mobile genetic elements and can be easily acquired through horizontal gene transfer. The ubiquity of different homologous and non-homologous TA systems within a single bacterial genome raises questions about their potential cross-interactions. Unspecific cross-talk between toxins and antitoxins of non-cognate modules may unbalance the ratio of the interacting partners and cause an increase in the free toxin level, which can be deleterious to the cell. Moreover, TA systems can be involved in broadly understood molecular networks as transcriptional regulators of other genes' expression or modulators of cellular mRNA stability. In nature, multiple copies of highly similar or identical TA systems are rather infrequent and probably represent a transition stage during evolution to complete insulation or decay of one of them. Nevertheless, several types of cross-interactions have been described in the literature to date. This implies a question of the possibility and consequences of the TA system cross-interactions, especially in the context of the practical application of the TA-based biotechnological and medical strategies, in which such TAs will be used outside their natural context, will be artificially introduced and induced in the new hosts. Thus, in this review, we discuss the prospective challenges of system cross-talks in the safety and effectiveness of TA system usage.

Structural insights of the toxin-antitoxin system VPA0770-VPA0769 in Vibrio parahaemolyticus.

Toxin-antitoxin (TA) systems are involved in both normal bacterial physiology and pathogenicity, including gene regulation, antibiotic resistance, and bacteria persistence under stressful environments. In pathogenic Vibrio parahaemolyticus, however, TA interaction and assembly remain largely unknown. In this work, we identified a new RES-Xre type II TA module, encoded by gene cluster vpa0770-vpa0769 on chromosome II of V. parahaemolyticus. Ectopic expression of the VPA0770 toxin rapidly arrests the growth of E. coli cells, which can be neutralized by co-expression of the VPA0769 antitoxin. To decipher the action mechanism, we determined the crystal structure of the VPA0770-VPA0769 TA complex. VPA0770 and VPA0769 proteins can assemble into two types of large complexes, a W-shaped hetero-hexamer and a donut-like hetero-dodecamer, in a concentration-dependent manner in solution. Disruption of the TA interface results in a loss of the antitoxic phenotype. The toxicity of the VPA0770 toxin, which harbors a NAD+-binding pocket, may be largely ascribed to its highly effective capability to degrade intracellular NAD+. Our study provides a structural basis for a better understanding of diverse molecular mechanisms employed by human pathogens.

A ParDE toxin-antitoxin system is responsible for the maintenance of the Yersinia virulence plasmid but not for type III secretion-associated growth inhibition.

Many Gram-negative pathogens utilize the type III secretion system (T3SS) to translocate virulence-promoting effector proteins into eukaryotic host cells. The activity of this system results in a severe reduction of bacterial growth and division, summarized as secretion-associated growth inhibition (SAGI). In Yersinia enterocolitica, the T3SS and related proteins are encoded on a virulence plasmid. We identified a ParDE-like toxin-antitoxin system on this virulence plasmid in genetic proximity to yopE, encoding a T3SS effector. Effectors are strongly upregulated upon activation of the T3SS, indicating a potential role of the ParDE system in the SAGI or maintenance of the virulence plasmid. Expression of the toxin ParE in trans resulted in reduced growth and elongated bacteria, highly reminiscent of the SAGI. Nevertheless, the activity of ParDE is not causal for the SAGI. T3SS activation did not influence ParDE activity; conversely, ParDE had no impact on T3SS assembly or activity itself. However, we found that ParDE ensures the presence of the T3SS across bacterial populations by reducing the loss of the virulence plasmid, especially under conditions relevant to infection. Despite this effect, a subset of bacteria lost the virulence plasmid and regained the ability to divide under secreting conditions, facilitating the possible emergence of T3SS-negative bacteria in late acute and persistent infections.

Investigating the role of Bacillus subtilis type II toxin-antitoxin system in drought stress survival.

Toxin-antitoxin (TA) systems, present in plasmids and bacterial chromosomes, are widespread in bacteria such as Bacillus subtilis and are known to be involved in growth regulation, bacterial tolerance to environmental stress conditions as well as biofilm formation. The aim of the current study was to investigate the role of TA systems in drought condition stress in B. subtilis isolates. The presence of TA systems including mazF/mazE and yobQ/yobR in B. subtilis (strain 168) was investigated using the polymerase chain reaction (PCR) method. TA system expression at 438 and 548 g/L of ethylene glycol concentrations was evaluated using real-time PCR method and sigB gene was used as internal control. The expression rate (fold change) of mazF toxin gene treated with 438 and 548 g/L of ethylene glycol was 6 and 8.4, respectively. This indicates an increase in the expression of this toxin in drought stress condition. Also, the fold change of mazE antitoxin in the treatment with 438 and 548 g/L of ethylene glycol was 8.6 and 5, respectively. While yobQ/yobR showed a decrease in expression in 438 and 548 g/L of ethylene glycol concentrations. So that the highest expression reduction (8.3) was observed for yobQ gene at the concentration of 548 g/L of ethylene glycol. Results of this study revealed the significant role of B. subtilis TA systems in drought stress which can be considered as the resistance mechanism of this bacterium under stress conditions.